First-principles calculations of spin and angle-resolved resonant photoemission spectra 

of Cr(llO) surfaces at the 2p - 3d resonance 
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A first principles approach for spin and angle resolved resonant photoemission is developed within 
multiple scattering theory and applied to a Cr(llO) surface at the 2p-3d resonance. The resonant 
photocurrent from this non ferromagnetic system is found to be strongly spin polarized by circularly 
polarized light, in agreement with experiments on antiferromagnetic and magnetically disordered 
systems. By comparing the antiferromagnetic and Pauli-paramagnetic phases of Cr, we explicitly 
show that the spin polarization of the photocurrent is independent of the existence of local magnetic 
moments, solving a long-standing debate on the origin of such polarization. New spin polarization 
effects are predicted for the paramagnetic phase even with unpolarized light, opening new directions 
for full mapping of spin interactions in macroscopically non magnetic or nanostructured systems. 

PACS numbers: 78.20.Bh,78.70.-g,75.20.Ls,79.60.-i 



In recent years, the theoretical description of absorp- 
tion/photocmission spectroscopy in the X-ray region has 
been boosted by the merge of density functional theory 
(DFT) with many body approaches such as dynamical 
mean field theory many body perturbation theory 

3h5( and by the development of time-dependent DFT 
6J. However, second order processes, like resonant inelas- 
tic X-ray scattering (RIXS) and resonant photoemission 
(RPES), remain a major challenge for theory. For RPES, 
existing approaches are semiempirical |7Hlfj|. based on a 
well defined two- holes final state and on small clusters, 
and thus do not take into account the derealization of 
intermediate states, the bandstructure of the system and 
multiple scattering effects in the propagation of photo- 
electrons. 

The huge experimental output from RPES on corre- 
lated materials @, IIlT - fl6| and the intriguing quest for a 
determination of local magnetic properties put forward 
by pioneering experiments [14|-|l6j call for advancements 
in the theoretical description of this spectroscopy. In ex- 
periments on CuO and Ni, it was shown that the RPES 
photocurrent with circular polarized light is spin polar- 
ized in antiferromagnets 14, [To| and Curie paramag- 
nets 161. It was claimed that a specific combination of 
spin resolved spectra provides a direct measure of the 
local magnetic moments 14Hl6l|. The issue is of funda- 
mental importance in the search for a tool to access the 
local magnetic properties in antiferromagnetic, magneti- 
cally disordered and/or nanostructured systems at their 
crossover with the transition temperature. The inter- 
pretation was however rejected on the basis of symmetry 
analysis 17J , but explicit calculations predicting the line- 
shape and intensity of such fundamental signal are still 
lacking and remain highly desirable. 

In this letter, we present the first ab-initio method 
for RPES in solids, based on a combined formulation 
within the real space multiple scattering (RSMS) ap- 



proach ill lj| and DFT, and its application to Cr(llO) 
at the 2p-3d resonance. By comparing the antiferromag- 
netic (AFM) and Pauli-paramagnetic (PM) phase of Cr, 
we solve the long-standing debate about the possibility 
to determine local magnetic moments in macroscopically 
non magnetic systems by means of spin resolved RPES 
with circular polarized light. New interesting effects in 
the PM phase by unpolarized light suggest that other 
mechanisms are active and could be exploited for map- 
ping the origin of the different spin polarization (SP) 
components in paramagnets and magnetically disordered 
systems. 

Theoretical formulation. The cross section for valence 
band photoemission to a final state \v, k) , where v de- 
notes a valence band hole and k a photoelectron state, is 
given by 



I(uj,q,k) = ^2 \T kv {^,q)\ 2 S(e k 



tvui) 



where fko and q are the photon energy and polariza- 
tion. Here the independent particle approximation has 
been assumed (i.e., all many-electron eigenstates are 
single Slater determinants corresponding to the same 
effective one-electron hamiltonian) . According to the 
Heisenberg-Kramcrs formula [20], the transition matrix 
element Tk v (uj,q) is the sum of a direct and a resonant 
term. In the latter, photon absorption leads to an inter- 
mediate state |c, u), with a core hole (c) and an electron in 
a formerly unoccupied state \u), which decays to the final 
state \v, k) through a participator Auger process 2ol 21|. 
To lowest order in the autoionization process, the transi- 
tion matrix element is given by 



Tkv(u,q) = (k\D q \v) + 



(kc\V(\vu) — \uv}) 
fvjj + e r — e„ — iT 



u\D q \c) 



(1) 

where D q is the dipole operator, V the Coulomb oper- 
ator and r the width of the intermediate state. Spec- 
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tator Auger decay leads to different, namely two-hole fi- 
nal states and is not considered here. Participator and 
spectator channels can in principle be separated experi- 
mentally by using a photon bandwidth smaller than the 
core-hole lifetime, as they show different photon energy 
dependence (linear for the participator, and no photon 
energy dependence for the spectator). Here we focus on 
the physical effects at the origin of spin polarization and 
dichroism as well as their directional-dependence in the 
"pure" participator channel. 
The RPES intensity can be written in a compact form 

as 



l(u),q,k) 



ijLL 1 a 
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Here, i, j label atomic sites, L = (Im) angular momentum 
and a spin quantum numbers. e v = e k — tko is the energy 



of the valence hole. The quantity IV L , 
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T ILL' 

is the essentially imaginary part of the scattering path 
operator. It comes from the simplification of the sum over 
delocalized valence states through the so called optical 
theorem in RSMS [22| and it contains the bandstructure 
information. The matrix elements M^£^(k) are given by 



■ {k)Aj L ^ iL (e k (Jk, f-vo) 



The Bjli (k) are the key quantities in the RSMS approach 
and represent the multiple scattering amplitudes of the 
continuum state k = (kcr^) [22| . The matrix elements 
AjL' ,iL(£k&k, £v&) are given by the sum of the direct ra- 
diative process (A D ), the resonant process with direct 
Coulomb decay (A c ) and the resonant process with the 
exchange decay (A x ), see Eq. ([T|). A D and A c are site- 
and spin-diagonal (~ 6ij6 ak a)- We have 



A D = (ie k L'<r\D q \ie v L<j) 
A c = - V / ,/,„ 



j' cL u L' u <7 u 
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The sums over unoccupied states u have been again 
simplified through the optical theorem. The exchange 
term A x is not strictly site-diagonal because of the non- 
locality of the exchange interaction together with the de- 
localized nature of the states u. In the RSMS approach 
the Coulomb matrix elements (fcc|T^|wu) and (kc\V\uv) 
can be exactly developed in one- and two-center terms. In 
metallic Cr, the Coulomb interaction is strongly screened. 
As a result, two-center terms are by at least one order 
of magnitude smaller than the one-center terms [23[ and 



have been neglected here. In general, the 2p-3d excited 
intermediate states might display excitonic effects, which 
could be taken account for with a Bethe-Salpeter descrip- 
tion 0, For Cr metal, these effects are quite small 
because of the large 3d band width (~ 7 eV) and effi- 
cient metallic screening of the core hole by nearly free 
Asp electrons, and thus neglected here. 

Photoemission spectra from Cr(llO) are calculated in 
RSMS with a cluster of 151 atoms (see Fig. la) and self- 
consistent spin polarized potentials, obtained by a scalar 
relativistic LMTO [24[ calculation for bulk Cr in the lo- 
cal spin density approximation. Except for the 2p core 
level, all states entering the RPES calculation are de- 
veloped in RSMS. The 2p orbital is obtained by solv- 
ing the scalar relativistic Schrodinger equation with self- 
consistent spin-polarized LMTO potentials. The 2p 3 / 2 
spin-orbit coupled states are then constructed using stan- 
dard angular momentum algebra and the spin-orbit cou- 
pling constant is taken from an atomic calculation [25j |. 
We consider the AFM order of CsCl-type which is a 
good approximation to the true spin density wave (SDW) 
ground state of Cr. The calculated magnetic moment 
is 0.74 fiB in reasonable agreement with experiment 
(0.62 (j,b)- At the (110) surface, the transverse SDW 
propagates along [100] or [010] [26[. Therefore, we take 
e z = [001] as magnetization and spin-quantization axis 
throughout this paper. We also consider the Pauli PM 
state, corresponding to a non-magnetic calculation. Spin 
orbit (SO) coupling of the valence and continuum states 
is neglected (it is as small as 0.03 eV for Cr-3d (27|V 

Results. The electronic structure of Cr(110) is well ac- 
counted for in the RSMS approach as can be seen from 
the comparison between the local density of states (DOS) 
of a Cr atom in the cluster and of bulk Cr (Fig.QJi). Non- 
resonant angle-resolved photoemission spectra (ARPES) 
are shown in Fig. QJ. Differences with respect to exper- 
iments [28| are expected as our approach does not con- 
tain local many-body interactions and layer-dependent 
potentials, which could play a role for a quantitative de- 
scription of the peak renormalization and dispersion be- 
haviour of the energetic structures [29]. However, the 
main features of the experimental spectra are reproduced 
in the calculation, confirming that RSMS provides a rea- 
sonably good description of valence band photoemission 
from metals as previously shown for Cu(lll) [22| . 

Spin resolved, angle integrated PES and RPES spec- 
tra are shown in Fig. [2] for the AFM phase and several 
photon energies across the L3-edge absorption threshold. 
Left circular polarized light incident along the magne- 
tization axis [001] is considered. In this "parallel" ge- 
ometry the spectra, right polarized light produces the 
same spectra but with up and down spin exchanged. The 
maximum peak intensity as a function of photon energy 
is plotted in Fig. [2Jd and shows the expected Fano pro- 
file. The first photon energy (551.0 eV) is too low to ex- 
cite the core electron and so only direct PES is possible. 
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FIG. 1: (a) Cr(110) cluster used in the RSMS calculations. 
The two magnetic sublattices of the AFM state are in red and 
blue, (b) DOS in the AFM phase for a bulk atom (LMTO) 
and a central atom in the cluster (RSMS). (c) ARPES spec- 
tra from Cr(110) along the (001) azimuth for different polar 
angles 9 with respect to the surface normal. Unpolarized light 
along the [001] axis was considered. Experimental data from 
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FIG. 2: a) Spin-resolved, angle integrated RPES and PES 
spectra of AFM Cr(110) with circular polarized light incom- 
ing along the spin quantization axis [001] and photon ener- 
gies across the Z/3-edge resonance. A gaussian broadening of 
0.27 eV FWHM was applied. Note the different intensity scale 
for hu = 554.4 eV. In PES, spin-up and down intensities are 
equal in all cases, b) Maximum peak intensity as a function 
of photon energy. 



When the photon energy is raised to 552.4 eV, just below 
the absorption edge, direct and resonant process interfere 



destructively, giving rise to the dip in the Fano profile. 
Strong resonant enhancement is observed between 552 
and 554.5 eV (see e.g. the spectrum for 554.4 eV), which 
corresponds to transitions from the 2p 3 / 2 level into the 
unoccupied Cr 3d band. At hv = 585.1 eV, well above 
threshold, the resonant spectrum goes back to the non- 
resonant one. 

The direct PES signal is non spin-polarized as expected 
for the AFM phase. Appreciable spin-polarization is, 
however, found in RPES. This effect is here obtained for 
the first time through first-principles calculations, and 
confirms the experimental finding in CuO 14j, that in 
AFM systems RPES at the 2p 3 / 2 -3d resonance is spin- 
polarized when circular polarized light is used. 

We now turn to angle and spin resolved spectra at 
maximum resonance (hj/=554.4 eV), focusing on their 
four "fundamental" combinations (and their relation 
to local magnetic properties), constructed by different 
choices of photoelectron spin (t,4-) and light helicity 
(+, -)=(left,right): 



tot = (t +) + (t -) + a 

spr = (t +) + (t -) - a 

die = (t +) - (t -) + a 

mix = (t +) - (t -) - a 



(i -) (total) 
(4- — ) (spin-resolved) 
(i — ) (dichroic) 
(| — ) (mixed) 



The "mixed" spectrum was the one considered in 
Refs 14 , [lfl and claimed to be sensitive to local mag- 
netic moments in non-ferromagnetic samples. 

The normal emission RPES spectra (Fig 3a, b) (total 
spectra) for parallel geometry consist of a single peak at 
0.8-0.9 eV binding energy, very similar to the low energy 
non resonant spectrum in Fig.lc (9 — 0°). AFM and PM 
spectra are almost identical except for a small shift of ~ 
0.1 eV, which reflects the small exchange splitting of the 
AFM Cr-3d bands. The dichroic (die) and spin-resolved 
(spr) signals vanish for both PM and AFM phase, as 
expected since the system is globally non-magnetic in 
both cases, and the set up is non chiral. 

However, the mixed signal is non-zero with a large am- 
plitude (~ 1/3 of total), in agreement with the experi- 
mental results in AFM CuO [14| . Surprisingly, we find a 
non-zero mixed signal not only in the AFM, but also in 
the PM phase with nearly the same intensity. It is impor- 
tant to note that we are not considering a Curie param- 
agnet (such as Ni above Tq [ill) with disordered and/or 
fluctuating magnetic moments, but a Pauli PM state, 
where the magnetization is strictly zero in all points of 
space. Therefore, our finding that the mixed signal is 
essentially unchanged when going from the AFM to the 
PM state unambiguously proves that it is unrelated to lo- 
cal magnetic moments, in contrast to the interpretation 
in Refs 0, 0. 

Rather than being of magnetic origin, the non-zero 
mixed signal is in fact induced by angular momentum 
transfer from the light helicity to the electron spin via 
SO in the core shell together with a strong exchange 
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FIG. 3: Angle-resolved fundamental spectra of Cr(llO). 
RPES as thick lines for hv = 544.4 eV. Normal (i.e., non 
resonant) ARPES as thin lines, intensity x 1000. All spec- 
tra are rescaled to equal peak height of RPES-tot. (a) 
AFM, (b-d) PM phase, (a-c) Normal emission, (d) Emis- 
sion in the xy-plane, off-normal by 23° (vector e in Fig.l a). 
Light incidence parallel (a,b,d) or perpendicular (c) to spin- 
quantization axis e 2 . Light vector in (c) is shown as p in 
Fig. la. 



effect in the decay process. To see this, consider light 
with left (+) helicity and a non-magnetic ground state. 
The 2p 3 / 2 -3<i optical transition has a larger amplitude for 
spin-up than for spin down electrons because of the dom- 
inantly parallel alignment of spin and orbit in 2p 3 / 2 - For 
example, for an empty or spherically symmetric 3d shell 
the intensity ratio is 5:3. Consider now a spin- up elec- 
tron transition. The RPES intermediate state has one 
extra spin-up electron in the 3d-shell (denoted uf) and a 
2p-holc of dominant spin-up character. This state decays 
through Coulomb interaction to the photoemission final 
state with one 3cZ-hole and the photoelectron. The direct 
Coulomb matrix elements is of the form (kcr, c\\V \vcr, u'l) 
which is independent of the photoelectron spin a. So the 
direct decay alone would lead to a spin-balanced pho- 
tocurrent. For the exchange decay, the matrix element 
is (kcr, cf\V\uf, vcr) ~ <5(er,t)- This is roughly as a large 
as the direct Coulomb term for spin-up electrons (the ra- 
dial matrix elements are exactly the same) but it is zero 
for spin-down electrons. Since the exchange matrix ele- 
ments are substracted from the direct terms in Eq. (1), 
the transition probability for spin-up electron emission is 
strongly reduced by the exchange process. This shows 
that a core- valence transition of a spin up electron leads, 
through autoionization, to a strongly spin polarized pho- 
tocurrent with a majority of spin down electrons. As 
mentioned before, left circular polarized light promotes 
dominantly spin- up electrons in the 2p 3 / 2 -3d transition. 



Therefore it produces a majority of spin-down photoelec- 
trons. Under the assumption of complete cancellation be- 
tween direct Coulomb and exchange matrix elements for 
parallel spins and by neglecting the direct valence photoe- 
mission, the ratio of spin-down to spin-up photoelectrons 
is 5:3, which corresponds to a spin-polarization (ratio of 
mixed over total signal) of —1/4. In angle integrated 
RPES at maximum resonance (Fig. 2a, /iz/=554.4 eV) we 
find a SP of —0.21, in good agreement with such model 
estimation. These values agree also well with the mea- 
sured spin-polarization in CuO [l4j and Ni [ltij . which 
is 10-40% depending on binding energy. Our findings 
clarify the physical mechanism inducing the presence of 
the mixed signal in both phases, and point to a critical 
re-examination of experimental observations. 

Interestingly, we find that, contrary to the previous 
set up, it is possible to have a net spin polarization sig- 
nal on the PM phase. This is possible under appropriate 
geometrical conditions, and even with unpolarized light. 
Such SP can be of opposite sign and be due to different 
active mechanisms. In Fig. 3c, normal emission spectra 
are shown for light incident along [110], i.e. perpendic- 
ular to the spin-quantization axis s=e z (perpendicular 
geometry). As before, the dichroic signal is zero, as light 
incidence (p) and electron emission vector (n) lie in a 
mirror plane of the surface (see Fig. 1 a). However, the 
set up (including spin resolution) is chiral, since the three 
vectors p, n and s form a right-handed frame. Thus 
SO-induced SP cannot be ruled out by symmetry and a 
small, positive SP (in this case transverse to the scatter- 
ing plane) is indeed observed in RPES, even for unpolar- 
ized light. A similar SP from PM surfaces for unpolar- 
ized light was theoretically predicted in direct PES [13] 
in a relativistic approach and confirmed by experiments 
was ascribed to broken symmetry due to the 



off-normal light incidence together with SO in the ini- 
tial states and phase shift differences. We do not observe 
this effect in non-resonant PES since the SO coupling in 
the Cr 3ci valence states is very weak and neglected here. 
However, for RPES, such SP has to be related to the dy- 
namical SP studied in atomic physics, which is known to 
be related to phase shift differences in the final outgoing 
waves, and to be generally small [HI, |34|. Our result con- 
firms that such SP exists for an atom embedded in a solid 
and that it survives to the multiple scattering effects. 

A SP signal in the PM phase is also present for paral- 
lel geometry with off- normal emission (Fig. 3d). In this 
case, the system composed by the surface, light incidence 
(along e z ) and electron emission vector, is chiral. There- 
fore a dichroic signal is observed even in non-resonant 
PES, known as circular dichroism in angular distribu- 
tion (35j . In RPES, the angular momentum of the photon 
is partly transferred to the electron spin through the SO 
coupling in the 2p shell, leading to non-zero intensity also 
for spin resolved and mixed signals. The spin polariza- 
tion is negative, i.e. photoelectrons are mainly polarized 
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antiparallel to their emission direction, because of the ex- 
change process in the autoionization decay. This finding 
suggests a Fano-like effect in resonant processes for off 
normal emission directions, which could be well studied 
along the same lines as direct PES on paramagnets [36| . 

In conclusion, we have presented a first-principles ap- 
proach for RPES in solids and its application to Cr(llO). 
By comparing Pauli PM and AFM states, we have shown 
that the mixed signal is essentially independent of local 
magnetic properties and we have clarified its origin: con- 
trary to previous interpretations, this effect is induced by 
an angular momentum transfer from the photon to the 
electron spin, through SO coupling in the core level and 
the exchange process in the autoionization decay. Our re- 
sults show that caution must be taken in linking the spin 
polarized or mixed signal to local magnetic moments, all 
the more so as the photoelectron spin may have compo- 
nents along and across the light helicity. New effects in 
the SP suggest that a mapping of spin interactions in 
paramagnets and disordered magnetic structures could 
be obtained via full tomography experiments at the core 
resonances even with unpolarized light. 
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